Is nonperturbative inflatino production during preheating a real threat to cosmology? 
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We discuss toy models where supersymmetry is broken due to non- vanishing time- varying vacuum 
expectation value of the inflaton field during preheating. We discuss the production of inflatino the 
superpartner of inflaton due to vacuum fluctuations and then we argue that they do not survive 
until nucleosynthesis and decay along with the inflaton to produce a thermal bath after preheating. 
Thus the only relevant remnant is the helicity ±3/2 gravitinos which can genuinely cause problem 
to nucleosynthesis. 
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Inflation is perhaps one of the best paradigms of the 
early Universe which solves some of the nagging prob- 
lems of the standard Big Bang cosmology One of 
the consequences of inflation is that it leaves the Uni- 
verse extremely cold, virtually devoid of entropy. Thus, 
the Universe requires to be reheated to a temperature 
at least more than O(MeV) to keep the successes of 
the Big Bang nucleosynthesis. Perhaps, one can imag- 
ine that the Universe reheats via the decay of the scalar 
field whose potential has dominated the Universe during 
the inflationary regime. Inflation leaves the inflaton field 
extremely homogeneous except for the quantum fluctua- 
tions produced during inflation. The perturbations keep 
their imprint intact to match the observed anisotropy in 
the present Universe which is one part in 10 5 ||. Once 
inflation ends, the mass of the inflaton field dominates 
over the Hubble rate of expansion, and the homogeneous 
inflaton field oscillates coherently around the bottom of 
the potential. If we assume chaotic inflation with a mas- 
sive inflaton field, to, and potential V = m 2 (j> 2 /2, then, 
during the coherent oscillations the average pressure of 
the Universe within one Hubble time vanishes over many 
oscillations. As a result the decaying energy density of 
the Universe behaves as if it were in a matter dominated 
era with = <p 2 /2 + m 2 (j> 2 /2 ~ a~ 3 , where a is the scale 
factor of the expanding Universe. After couple of oscilla- 
tions the energy density in the scalar field redshifts away 
in the same way as in the pressureless fluid but this does 
not lead to a radiation dominated Universe. To obtain a 
radiation dominated era, the inflaton field has to decay to 
other particles which will eventually lead to a thermalized 
plasma with a finite temperature, usually known as the 
reheat temperature of the Universe. The inflaton decays 
when the Hubble parameter H ~ T^, where is the de- 
cay rate. The decay rate essentially depends on the kind 
of couplings the inflaton has to other particles |3) . How- 
ever, in between the end of inflation and the beginning 
of the radiation era there can be an explosive produc- 
tion of particles purely due to non-thermal effects. This 
new wisdom has been realized in Rcfs. 01. This is due 



to the fact that the oscillations in the inflaton field are 
extremely coherent and act as a Bose condensate fluid. 
So, in principle, one can study the quantum fluctuations 
of the inflaton quanta as well as bosonic and fermionic 
fields which are coupled to the inflaton field via Yukawa, 
gauge, or non-renormalizable couplings. Effectively, the 
problem turns out to be quantizing the bosonic and the 
fermionic fields in a time-varying inflaton background. 
This leads to an explosive production of particles which 
does not depend on the background temperature and it 
is purely an offshoot of a non-perturbative analysis. The 
production of bosons and fermions differs in its nature 
due to Pauli's exclusion principle, which prohibits ex- 
cessive production of fermions compared to their bosonic 
counterparts || . In this regard, recently it has been real- 
ized that like fermions with spin 1/2, other fermions with 
higher spin can also be created from the vacuum fluctu- 
ations in a time- varying scalar background. In Ref. 0, 
the authors have noticed that inspite of Planck mass sup- 
pressed couplings of spin 3/2 particles to other fields, it is 
possible to excite them due to vacuum fluctuations. This 
has lead to many consequences which we briefly discuss 
in the next paragraph. 

The spin 3/2 gravitino occurs in supersymmetric theo- 
ries as a superpartner of the graviton. A massive spin 3/2 
has four helicity states ±3/2 and ±1/2. A massless grav- 
itinos only possess ±3/2 helicity states. However, once 
supersymmetry is broken the gravitinos become massive, 
and they possess all four helicity states. In the early 
Universe supersymmetry can be broken due to non-zero 
vacuum contribution of the inflaton energy density. If the 
inflaton field is a scalar component of a chiral multiplet, 
then spontaneous supersymmetry breaking due to F-term 
leads to non-zero expectation value of the fermionic field 



(0|^|0) = (-i m-e G/2 G*0 ^0. 



(1) 



where £ is the infinitesimal Grassmann-odd parameter, <j) 
is a scalar field responsible for inflation, <f> is the fermionic 
component of the inflaton in a single chiral field model, 
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which we may call inflatino, G is the Kahler function 
defined below, and is the derivative of the Kahler 
function with respect to the inflaton field. The Kahler 
function is defined by 

<k<t f\w\\ f . 

G= ^- +U [w) ' (2) 

where we have assumed minimal Kahler function, W is 
the superpotential and M = M p /\f8n is the reduced 
Planck mass. Out of the two terms present in the right- 
hand side of Eq. ([!]), we notice that the first term gives a 
non-zero contribution during and after inflation, particu- 
larly during preheating. Therefore, the dynamical effects 
of the inflaton field breaks supersymmetry. 

Soon after it has been realized that the helicity ±3/2 
states of a massive gravitino can be produced non- 
perturbatively |6), it has been shown that the helicity 
±1/2 states of a massive gravitino can also be produced 
from vacuum fluctuations B. However, they are more 
abundantly produced compared to that of the helicity 
±3/2 states. This can be easily understood in a simple 
way. For the creation of particles from vacuum fluctua- 
tions, the adiabaticity condition has to be broken which 
is usually measured by a rate of change of a time- varying 
frequency of a given momentum mode. For fermions the 
frequency depends on an effective mass parameter. For 
example, for the helicity ±3/2, the mass parameter is es- 
sentially Planck mass suppressed, ft has been noticed in 
Refs. that the helicity ±1/2 states are massive due to 
the fact that they eat the mass of the fermionic compo- 
nent of the inflaton. This statement is true for a single 
chiral case and it has been pointed out that for helic- 
ity ±1/2 gravitinos, the adiabaticity condition is broken 
much more strongly compared to that of helicity ±3/2 
gravitinos |p|-|l^|. 

All these results obtained were interesting because 
gravitino plays a key role in a standard Big Bang cos- 
mology. If supersymmetry is required to solve the 
gauge hierarchy problem, then, in the gravity medi- 
ated supersymmetry breaking, the gravitino gets a mass 
around, O(TeV). Since their couplings to other parti- 
cles are Planck mass suppressed, the life time of the 
gravitino at rest is quite long, r 3 / 2 ~ -^p/ m 3/2 ~ 
10 5 (m 3 / 2 /TeU)~ 3 sec, We know that successful nu- 
cleosynthesis depends on the baryon abundance: Yb(T < 
MeV) = n 3/2 /n 7 = 1CT 10 @. The gravitino decay 
products can easily change this ratio. Their decay prod- 
ucts such as gauge bosons and its gaugino partners, or 
high energy photons, can generate a large entropy which 
will heat up the photons compared to r and \x neutrinos. 
The abundance of neutrinos essentially determines the 
4 He abundance. It was first pointed out in Ref. Jl5j that 
the gravitino mass must be larger than ~ 10 TeV in order 
to keep the successes of the Big Bang nucleosynthesis. On 
contrary, if the gravitinos were stable, and if their mass 
exceeded 1 KeV, they could easily overdose the Universe 
in absence of inflation m&. However, after the end of in- 



flation the gravitinos can be produced from the thermal 
bath and this constraints the temperature of the thermal 
bath in order not to over produce them. At the time of 
nucleosynthesis the abundance is given in terms of the 
reheat temperature, Y 3/2 (T < MeV) - 10" 2 (T rh /Af p ) 
|l7| . Thus, we see that there exists a strong constraint 
on the reheat temperature, T r h < 10 10 GeV, in order to 
maintain the baryon abundance one part in 10 10 during 
nucleosynthesis. Since we know that non-perturbative 
creation of particles does not depend on temperature, it 
will be difficult to constrain a general parameter other 
than the model parameters. Hence, this leads to a nat- 
ural suspicion that perhaps non-perturbative production 
of helicity ±1/2 gravitinos will cause a problem to nucle- 
osynthesis. 

The important point is that the inflaton has to com- 
pletely decay to give rise to a thermal bath with a re- 
heat temperature at least more than a O(MeV), and the 
fermionic component of the inflaton, known as inflatino, 
inevitably decays along with the inflaton. We know dur- 
ing the inflaton oscillations, the helicity ±1/2 states of 
the gravitino eat the mass of the inflatino, and they es- 
sentially behave as an inflatino when the amplitude of 
the inflaton oscillations has considerably dropped below 
Mp. As a result they must also decay along with the in- 
flaton. As we shall see that this argument is quite robust 
and it should not depend if there were any other source 
of supersymmetry breaking other than the inflaton sec- 
tor. During the preheating era of the Universe it is quite 
natural to think that the supersymmetry breaking due 
to the energy denisty stored in the inflaton oscillations is 
far the most dominant source. 

We will begin with an introduction of a supersymmet- 
ric inflationary model with a single multiplet, and then 
we discuss decay rates of the inflaton and the inflatino in 
two models: namely with Planck mass suppressed cou- 
pling, and with Yukawa couplings to the visible sector. 
We then establish an equivalence between the helicity 
±1/2 gravitino interactions to its supercurrent to that of 
the inflatino interactions in the supergravity Lagrangian 
when the amplitude of the inflaton oscillations are small 
compared to the Planck mass. In the last section we give 
a qualitative discussion upon the gravitino decay when 
more than one chiral fields are present. 

I. MODELS WITH A SINGLE MULTIPLET 

In the most part of this paper we shall focus on mod- 
els where supersymmtery is broken by a single multiplet 
and also responsible for producing inflation in the early 
Universe. Nevertheless, to solve the low-energy (i.e. elec- 
troweak scale) supersymmetry breaking we may require 
some other sector, which can be a hidden sector, which 
we shall not take into account here. In our case the source 
of a time- varying supersymmetry breaking is the oscilla- 
tions in the inflaton field <f>. During these oscillations the 
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fermionic partner of the infiaton which we call here in- 
flatino, whose mass is equal to the mass of the infiaton, is 
eaten by the helicity ±1/2 components of the gravitino 
to produce a massive gravitino. It has been suggested 
by many authors in Refs. |^,|l^,^2), that in a limit when 
\<f>\ < M p , it is possible to use the inflatino mode equation 
to study the behavior of the helicity ±1/2 states of the 
gravitino. This can as well be understood from the point 
of view of an Equivalence Theorem (ET), which demands 
that when the energy scale E 3> m3/2(i), the wavefunc- 
tion of the gravitino for the helicity ±1/2 components is 
approximately proportional to p^/m^^it), where p^ is 
the momentum of the gravitino and m^i^t) is its time- 
varying effective mass. However, there is a word of cau- 
tion regarding the validity of ET in our calculation during 
the oscillations of the infiaton. In principle, the time- 
varying mass of the gravitino can be larger than the mo- 
mentum during the oscillations, or perhaps m^/2{t) oc p^, 
and in both the cases ET cannot be trusted during the in- 
fiaton oscillations. However, studying the inflatino mode 
equation is not futile, because when the amplitude of 
the oscillations die down due to the expansion of the 
Universe, it is possible to identify the high momentum 
Fourier modes of the inflatino with those of the helic- 
ity ±1/2 gravitino. Therefore, only in those regions we 
can identify the inflatino to the helicity ±1/2 gravitinos 
(in a Fourier space), and we can therefore identify the 
Bogolyubov coefficients which are related to the number 
density of the produced helicity ±1/2 gravitinos. In this 
paper we are going to argue that ET can also be used to 
study the decay of the helicity ±1/2 gravitinos. However, 
this means that by using ET we shall be able to match 
the coupling strength of the helicity ±1/2 gravitinos to 
that of the inflatinos. This we shall discuss in the coming 
sections. 



A. Infiaton decaying via gravitational coupling 

As a first example we consider a new inflation model 
proposed in Ref . [1£| . In this model the two distinct sec- 
tors are the infiaton sector and the visible sector. These 
sectors interact with each other only gravitationally, and 
can be considered separately in the superpotential. The 
construction of the infiaton sector demands supersym- 
metry is restored in the global minimum. While setting 
the cosmological constant to zero, the simplest form of a 
superpotential emerges Q 



I 



A 2 
M 



($-M) 5 



(3) 



where A determines the scale of inflation. Here we have 
denoted $ as a superfield in the infiaton sector. The 
amplitude of the density perturbations produced during 
inflation by the infiaton, 0, is fixed by the COBE scale, 
which constraints A/M as 5 x 1CU 3 . With this choice 
of superpotential, inflation occurs for <fi <§; M; the os- 
cillations take place around the minimum of potential 



4>o = M, with a frequency ~ A 2 /M. The scalar po- 
tential derived from the above superpotential has a form 



K = e £,(i*,-i/ M ) 2 



dW ti 



tW u 



M 2 



\Wu 



M 2 



(4) 



where we have assumed minimal Kahler function and we 
consider the total superpotential to be 



I 



(5) 



where L can be recognized as a visible sector which con- 
tains the light degrees of freedom. Before we begin our 
discussion on decaying infiaton, we mention some of the 
essential points related to this model. The dominant cou- 
pling of the infiaton to other light degrees of freedom can 
be read from the potential, Eq. Just by expanding 
the interference term in Eq. (^), we notice that the in- 
fiaton field can decay only via trilinear coupling to the 
scalars. This certainly prevents creation of such scalar 
fields via parametric resonance. Hence, the decay of the 
infiaton is essentially perturbative in nature. Under the 
condition > H , the decay rate of the infiaton does not 
depend upon the curvature of the universe. However, the 
infiaton field has a time varying amplitude, thus it is nec- 
essary to virialize the mean value of the field. Otherwise, 
we may expand the infiaton field around its minimum 
value M, by assuming 



4>' = (j) - M- , 



(6) 



where <fi is assumed to have pure oscillatory part with an 
amplitude much less than one in units of reduced Planck 
mass. With the help of Eq. (||), it is easy to evaluate the 
interference terms coming from the first squared term in 
the bracket in Eq. (||). The leading order term in the ex- 
pansion generates trilinear coupling to the matter sector 
from L with a gravitational strength ~ A 2 /M 2 , corre- 
sponding to a decay width ~ (A 2 /M 2 ) 2 . Since 
the mass of the infiaton is to^ ~ A 2 /M, this gives a fi- 
nite decay width of the infiaton p8| 



M 5 



(7) 



If we assume that the infiaton energy is converted into 
radiation according to 



P4> ~ 2~q3* T * > 



(8) 



where is the relativistic degrees of freedom, then the 
reheat temperature of the Universe can be estimated by 



(9) 
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For A/M ~ 5 x 10 3 , the reheat temperature is around 
T r - 10 s GeV. 

With this introduction we may now turn our attention 
to the decay of the hclicity ±1/2 gravitinos which are cre- 
ated during the oscillations of the inflaton field from the 
vacuum fluctuations. We remind that gravitino produc- 
tion is completely non-thermal, and we cannot associate 
their number density to any particular thermal bath. We 
also notice that the mass of the gravitino need not nec- 
essarily be that of a gravitino mass around O(TeV). Es- 
pecially, if the inflaton sector has a supersymmetric pre- 
serving minimum with a zero cosmological constant, then 
the mass of the gravitino vanishes after the end of reheat- 
ing, provided there is no other source of supersymmetry 
breaking in any other sector. However, the situation can 
be little bit different if there exists other sources of super- 
symmetry breaking. This we shall briefly discuss in the 
next section. In the present section we shall concentrate 
upon the former case where we argue that whatever he- 
licity ±1/2 gravitinos are created during preheating shall 
have to decay along with the inflaton to have a successful 
Big Bang nucleosynthesis []. 

The equation of motion for the helicity ±1/2 gravitino 
in a cosmological background has been derived in the lit- 
erature by using alternative approaches |^-[l^]. The im- 
portant realization is that when the amplitude of the os- 
cillations is much smaller than the reduced Planck mass, 
the equation of motion for the hclicity ±1/2 gravitino 
is effectively that of the goldstino in a global supersym- 
metric limit. For a single chiral case the goldstino is 
equivalent to the inflatino up to a phase. Here we sim- 
ply quote some of the useful formulae which have already 
been established in Refs. Pp|Jl2]|. The evolution of the 
inflatino, which we define here as </>, is given by || 

ij°4> - k<t> - m ctf 4> = , (10) 

where k = "/ l ki, and fcj are components of the physi- 
cal momentum, while 7* are the gamma matrices. The 



* Here we would like to point out that eventhough we are 
considering inflaton to be the only source of supersymmetry 
breaking which might seems unrealistic at first point but the 
analysis is much simpler in this case and our arguments hold 
true even if there is another source of supersymmetry break- 
ing. It is possible to consider a Polnyi sector which breaks 
supersymmtry in the hidden sector of the theory but as we 
shall see in the next section that its mass contribution to the 
mass of the goldstino is very small and of the order of TeV. 
As a result the goldstino mass is essentially dominated by the 
inflatino mass. In such a circumstance our present analysis 
of single chiral field is quite general and as long as supersym- 
metry breaking due to the inflaton sector dominates over any 
other sector which is quite normal to think, we can apply our 
results faily well. All that we require is that the helicity ±1/2 
gravitino mass is essentially being contributed by the inflatino 
mass. 



validity of the above equation holds only in a global su- 
persymmetric limit. 

When the amplitude of the inflaton oscillations \<j>\ <C 
M, the effective mass of the helicity ±1/2 gravitinos, for 
a single chiral field and after phase rotation of the helicity 
±1/2 gravitino field, is simply the mass of the fermionic 
component of the inflaton field, which yields Q 

d 2 l 

m cS = — , (11) 

where I is the inflaton superpotential. For a simple su- 
perpotential Eq. (^), the effective mass for the helicity 
±1/2 gravitino turns out to be equal to 

A 2 

m±i/ 2 ^ TOefi ~ — , (12) 

which is the same as the mass of the inflaton |L8| . On the 
other hand, for the same superpotential the other helicity 
±3/2 gravitinos have mass given by 

™±3/2 = e yp m{1a) ' (13) 

where we have assumed that the visible sector L does not 
contribute to the gravitino mass. This is quite apparent 
from the above expressions Eqs. ([l^ , |l3| ), that the mass 
of the helicity ±3/2 gravitinos is not only suppressed by 
the reduced Planck mass, but it also contains time vary- 
ing amplitude of the oscillations; ~ 4>(t), which becomes 
vanishingly small near the bottom of the potential. This 
is quite obvious because mass of the helicity ±3/2 grav- 
itinos is essentially generated by the dynamics of the in- 
flaton field, and it must vanish when supersymmetry is 
restored at the global minimum of the potential. Before 
we begin our discussion on the decay of gravitinos, we 
compare different mass scales with the Hubble expansion. 
For the superpotential Eq. (||), the Hubble parameter it 
is given by H ~ (A 2 /M)(0(i)/M). This leads to a simple 
inequality in various mass scales which we must bear in 
mind 

TO0 ~ TO±l/2 > H > TO± 3 /2 ■ (14) 



B. Inflatino interactions 

In this subsection we analyze the decay rate of the 
inflatino. We consider a following interaction which can 
be found in the matter Lagrangian H 

|det e\~ l C = -^eW&CPxiXfl. + h-c. , (15) 

where G % is the derivative of the Kahler potential with 
respect to left and right chiral components. We can fix 
the index; i = <fi, corresponding to the inflaton sector. 
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This leaves the other index j to run on the chiral com- 
ponents of the visible sector L. Now according to our 
previous discussion on the inflaton decay, here again, we 
may assume that the inflatino is decaying into particles 
and sparticles of the sector L with an interaction of the 
form (fcxj&ji where <j> is the inflatino, \j is the fermionic 
component, and 4>j is the bosonic fields belonging to the 
sector L. The decay is essentially mediated via gravita- 
tional coupling strength ~ A 2 /M 2 . This corresponds to 
a decay width of the inflatino with mass ~ A 2 /M, which 
yields 



M 5 



(16) 



This decay rate is the same as the decay rate of the in- 
flaton. The result is not surprising because the inflatino 
mass is similar to the inflaton mass, and both interact 
gravitationally to the visible sector particles. Now, if we 
argue that the helicity ±1/2 states of the gravitino essen- 
tially behave as an inflatino in a global supersymmetric 
limit, which corresponds to demanding that the ampli- 
tude of the inflaton oscillations \<j>\ <C M, then and only 
then, we may argue that the helicity ±1/2 gravitinos de- 
cay along with the inflaton into visible sector particles. 
This is an important and generic conclusion which bears 
more attention. Our statement is only true provided we 
believe in the equivalence between inflatino and the he- 
licity ±1/2 gravitinos at late stages of the inflaton oscil- 
lations, i.e. when Eq. (14|) is satisfied. 

Intuitively, our result makes sense, because if super- 
symmetry is restored at the bottom of the potential in 
the absolute minimum, then, only the ±3/2 components 
of the gravitino should survive, and not the ±1/2 com- 
ponents of the gravitino. Thus, the helicity ±1/2 states 
must decay along with the inflaton decay. This situation 
could have been different if there were a hidden sector, 
which was responsible for supersymmetry breaking at an 
intermediate scale, which would then communicated to 
the visible sector at low scale. This we shall discuss in the 
last section. So far we have studied only the inflatino in- 
teractions. However, to be more concrete we must study 
the gravitino interactions which we shall discuss in the 
next subsection. 



light degrees of freedom. We have neglected the vector 
multiplets in the above equation and assumed <fi to be 
homogeneous. The total derivative -D^is defined by 



D„ 



<9 M + -U^ab<J 



a b 



(18) 



where w^ab is the spin connection. 

It is to be mentioned that apart from the derivative 
coupling of the chiral field to the gravitino, we have an 
extra interaction term which is not usually considered 
otherwise. In fact the interaction terms proportional to 
7 A '\E' M are usually not necessary in a static limit of the 
background field (i.e. inflaton field ), because 7 Mx I , A1 = 
acts as a constraint for a gravitino field in a static back- 
ground. However, this need not be true in a non-static 
background. It has been shown that in an expanding Uni- 
verse, and in a time- varying inflaton background, ±1/2 
helicity states follow 7 M 1 4 r ' i ^ ||. Although, the same 
constraint continues to hold good for the helicity ±3/2 
components of the gravitino in the same background 
along with the Dirac equation Thus, both the terms 
in Eq. ( p"7j ) should be taken into account to study the 
efficient decay of the gravitino. In this subsection we will 
study the decay by assuming the validity of the equiv- 
alence between the helicity ±1/2 states of the gravitino 
and that of the inflatino at late stages of oscillations, i.e. 
when Eq. ( |l4| ) is satisfied. 

After several oscillations of the inflaton field \<f)\ <C M, 
or, equivalently H <C m. Note that under this condition, 
the kinematics of the inflaton, such as decay rate does 
not depend on the curvature of the Universe. As a re- 
sult of this the decay rate of the inflaton coicides with 
that of the flat space-time limit. All the fields whose 
effective mass is larger than the Hubble parameter dur- 
ing the oscillations of the inflaton would actually not feel 
any effect of curvature of the Universe. Since, the effec- 
tive mass of the helicity ±1/2 gravitino is similar to the 
mass of the inflaton, and, it is much larger than the Hub- 
ble parameter, suggests that we can study their evolution 
by neglecting the curvature of the Universe. Therefore, 
we replace ±1/2 helicity of the gravitino by an ansatz 



2 M 

3 7^ 



(19) 



C. Interactions of the gravitino 

The gravitino interaction terms appear from the cou- 
plings between the gravitino field and the supercurrcnt 



4>J 



+h.c. , (17) 



where fi stands for the space-time index, Xi is a fermionic 
field and <p l is a bosonic field. Here the subscripts i,j 
correspond to the visible sector L, which contains the 



where rj represents the spin 1/2 fermionic field, which we 
shall interpret as a goldstino instead of inflatino. At this 
moment this prescription seems to be unwarranted, but, 
we shall see that this choice of derivative wavefunction 
leads to the interactions of the helicity ±1/2 gravitino to 
that of the inflatino. This prescription has also been used 
in Ref. Q. The goldstino is however related to the in- 
flatino by a phase factor, and, it is expressed in Eq. (|2^). 
The above expression is exactly the wavefunction of the 
helicity ±1/2 gravitino in terms of goldstino in the limit 
when the energy scale of the gravitino is larger than its 
effective mass. This expression has been previously used 
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in Refs. fll9|-pl"|], where the authors have been studying 
the scattering processes of the helicity ±1/2 gravitino in 
a limiting case when the energy scale of the gravitino is 
larger than its mass in a flat space-time. Here, again we 
have a similar situation where the helicity ±1/2 gravitino 
does not feel the Hubble expansion, however, the only dif- 
ference is that now supcrsymmetry is broken due to the 
oscillating scalar field rather than the static vacuum con- 
tribution. This obviously makes the difference. Notice, 
that instead of the mass of the gravitino m 3 / 2 ~ ITeV, 

1/2 

now we have pj /M, this is precisly because of the rea- 
son we have mentioned above, here supersymmetry is 
effectively broken due to the presence of the finite energy 
contribution of the inflaton field. The oscillations of the 
inflaton field also ensures that one has to include the in- 
teractions which are also proportional to 7 /i, 3> M . Another 
fact that the equation of motion of helicity ±1/2 grav- 
itino is the same as that of the goldstino for H <C m, as 
indicated in Refs. further supports our approxima- 

tion. We reiterate that we shall always work in a limit 
0/M< 1. 

Substituting Eq. @, in Eq. @, we derive an effec- 
tive Lagrangian. Which yields 



1 + 75 



X + h.c. , (20) 



where <p denotes the bosonic field in a visible sector L, 
the spinor component is defined by x; its mass can be 
written as 



expression we have neglected the time derivative of the 
energy density. Eq. (|23| ) can be further simplified if we 
assume that the mass splitting between ip and \ is due 
to supersymmetry breaking by the inflaton oscillations. 
To simplify the situation we will be assuming that the 
visible sector must contain the quadratic terms in the 
superpotential. This can be written as 



1 9 2 L 2 



(24) 



where A(ip,x) denotes the superfield, and ... terms can 
contribute due to other possibilities in the superpoten- 
tial, which we shall not take into account here. Now we 
will explicitly show that if the inflaton sector and the vis- 
ible sector interacts gravitationally, then it is possible to 
derive an effective inflatino Lagrangian which will have a 
similar coupling to Eq. ( |I1f ) . To get the desired result we 
first need to know the mass splitting between the fields 
of the visible sector. To get the mass splitting we expand 
Eq. (Q), with the help of Eqs. (||) and (p4|), while con- 
sidering only the dominant terms in the potential which 
is due to the interference terms. We also take help of 
Eq. @ to obtain 



«4r 



,0 W 1 dl 



(25) 



and, similar expression holds for m^j — m x , except for 
the negative sign on the right hand side. Similarly, one 
may also obtain 



dm 



Of 



(26) 



>72M : 



: d 2 L 
d^ 



(21) 



We also notice that the goldstino can be expressed as in 
Refs. 



We will use m x quite often to compactify our notations. 
In fact the mass of the fermion need not be a constant and 
can have field dependency. We can simplify Eq. ( |20| ) if we 
follow the below mentioned identities for the Majorana 
spinors 



vi^i b x = Xl^l^V , 

V15X = X15V ■ 



(22) 



With the help of Eq. ( p2| ) , we can derive an effective La- 
grangian after some algebraic manipulation, which reads 



£eff = 



1 



1/2 



2 2 \ • ^ m x 

m V R ~ m X l ^R- 7 ?* _ 1 



dt 

dm 



h.c. + total derivative , 



dt 



(23) 



1/2 



U <? _ 4> 2 /2M 2 dl 

dt dd> 



(27) 



where inflatino is represented by <f>. Here we have ex- 
plicitly used the fact that the dynamics of the inflaton 
is breaking supersymmetry. In Eq. (p7[), we have only 
retained the leading order terms and neglected C(l/M 2 ) 
terms. Now, with the help of Eqs. ( f25| , |26] , p7| ), we simplify 
Eq. pi) 



£ cS = m x — lpR( f )x + im x —(p I( j>j 5 x + h.c. , (28) 

where we have used the identity 7°^ = 7 , and the fact 
that the following relation holds 



7 dt d6 V d<b %1 dt 



. 0' 



~P4>- (29) 



where m v n denotes the real part of the light bosonic field 
ipn residing in the sector L. While deriving the above 



Now it is interesting to note that Eq. (28), upto a leading 
order actually leads to a familiar form 



G 



c, 



cff 



dcf> dip 



h.c. 



(30) 



which is the inflatino coupling in Eq. (|T^). While deriv- 
ing the above expression we have assumed Eq. . This 
clearly indicates that at late time of the inflaton oscilla- 
tions when we recognize the helicity ±1/2 component 
of the gravitino as a goldstino, we essentially get simi- 
lar coupling to the visible sector as that of the inflatino. 
This is the most important equivalence which establishes 
the fact that, since, for any successful inflationary model 
inflaton has to decay, and, so does the inflatino, the he- 
licity ±1/2 component of the gravitino must also decay 
if the inflaton oscillations is the only viable source of su- 
persymmetry breaking at that time. Our result is strictly 
correct for a single chiral field responsible for supersym- 
metry breaking. A further generalization to multi-chiral 
field supersymmetry breaking is more involved and we 
leave this for our future investigation. 

Now, we move onto a toy model where the inflaton 
sector and the visible sectors are coupled via Yukawa 
couplings. We will establish similar result as we have 
already obtained in this section. 



D. Model with a Yukawa coupling to the Inflaton 

As a second example we consider a model with a fol- 
lowing superpotential 



W= -m<$> 2 + -h<PY, 2 , 
2 2 



(31) 



where $ contains the inflaton field, which is responsi- 
ble for the slow-roll inflation. However, now the inflaton 
field has an explicit Yukawa coupling to the matter sector 
given by the second term in Eq. ([n]). Such a coupling will 
enable inflaton to decay much more efficiently. Such a 
superpotential leads to interaction terms hmcfxra, h(f>aa, 
h(f>aa, where <fi is the inflaton field, ip is the inflatino, a is 
a light bosonic field, and its fermionic partner has been 
denoted by a. The estimated rate of the inflaton decaying 
to fermionic component a a is given by Y ^ ~ (h 2 /8ir)m. 

In general the Yukawa coupling between $ and £ mul- 
tiplets can also result in the oscillations along the a field. 
If the a field eventually decays into other products much 
before the oscillations in a commences, then, it can still 
be a viable model to imagine that supersymmetry is 
broken by the inflaton field only. But in general, this 
may lead to a more complicated situation where super- 
symmtery is broken by several multiplets. However, it 
is possible to prevent this provided we require that the 
0-induced mass to the a field is much smaller than the 
Hubble expansion, i.e. hip < H, which implies h < m/M. 
We note that this will also insure that a and er are not 
produced via parametric resonance. A viable choice of 
parameters which can lead to an inflationary paradigm 
for the <fi field in a quadratic potential are; m = 10 
GeV, and, a small Yukawa coupling h = 10~ 7 , which 



ensures that at late stages of the inflaton oscillations, 
cj)/M < 10~ 14 , the inflaton is decaying perturbatively. 
Following our previous discussion, again, we argue here 
that since the inflatino mass is same as that of the mass 
of the inflaton, and, if the helicity ±1/2 components of 
the gravitino is recognized as inflatino at late stages of 
the inflaton oscillations, then they must decay to er, or er 
via a Yukawa coupling. 

So far, we have been looking upon direct inflatino cou- 
pling to er and a. However, we may now repeat the same 
analysis as we have shown earlier that indeed the helicity 
±1/2 components of the gravitino has a similar coupling 
as that of the inflatino by using the equivalence theorem. 
The generalization is quite simple and we recognize that 



d 2 W 
<9£ 2 



h<I>. 



(32) 



where £ and <& are the superfields denoted in Eq. Q31 
The mass of the fermion in this model is given by 



h((f>). Following our earlier argument we can find 



out the leading order contribution to the mass splitting, 
which yields 



(33) 



where I = (l/2)m$ 2 . Similarly, one can also derive an 
expression for m^j— m|, which differs from the above by a 
negative sign. While deriving Eq. (|3|), we have neglected 
the Planck mass suppressed contributions which would 
anyhow be insignificant at late times. The analogue of 
Eq. ( ^6| ) can be expressed as 



drria 



(34) 



and, one can also derive an effective Lagrangian with the 
help of Eqs. (||), ©, © and (|J). Which yields 



£ cff ~ h(a R (f>a + iai4>j 5 a) + h.c. , 



(35) 



where inflatino is denoted by <p. After some calculation 
it can be shown that Eq. (^5|), actually leads to an ex- 
pression 



C. 



cff 



ha* 



(36) 



This reinsures our earlier claim that the equivalence be- 
tween the helicity ±1/2 gravitino and the goldstino is vi- 
able at late times of the inflaton oscillations. This equiv- 
alence is not only important for studying the production 
of the helicity ±1/2 components of the gravitino, but also 
describing the decay of the helicity ±1/2 gravitinos. 

So far, we have not spoken any word on the other he- 
licity states of the gravitino, namely ±3/2. The reason 
is it is extremely difficult to study their decay, precisely 
because the mass of the helicity ±3/2 is solely due to the 
dynamics of the inflaton field |q| . Their effective mass is 
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Planck mass suppressed, and also depends on the ampli- 
tude of the oscillations of the inflaton field. This leads 
to an obvious result that if there is no other source of 
supersymmctry breaking other than the inflaton oscilla- 
tions, then, the effective mass for the helicity ±3/2 com- 
ponent should vanish at the end of reheating. It is dif- 
ficult to make a precise calculation for the decay of the 
helicity ±3/2 gravitinos. However, we believe that their 
survival does not depend on the inflaton decay as they 
have no goldstino nature. Next, we discuss qualitatively 
what would happen if the hidden sector supersymmctry 
breaking is also taken into account. 



II. MODELS WITH SEVERAL MULTIPLETS 

Once we invoke more than one sectors, and treat them 
at equal level, the problem of gravitino production be- 
comes more complicated. This problem has been ad- 



dressed in Refs. |9||lJ,|12j to some extent, and yet lot to 
be understood in this direction. In this case it has been 
realized that the goldstino is a linear combination of all 
the fermions, and as a result, even if we use the goldstino- 
gravitino equivalence we cannot in general guarantee that 
a major contribution to the goldstino mass is coming 
from the fermionic component of the inflaton. There are 
some interesting cases where the multi field case can be 
expressed as a single field, such as supersymmetric hybrid 
inflation model where effectively the two fields behave as 
if there were a single degree of freedom jnj. In such a 
model it is possible to extract the goldstino mass, which 
is again of the order of the inflaton mass. One can then 
discuss the decay rate of the inflatino in this model , 
and the inflatino decay rate to the light degrees of free- 
dom would exactly be the same as that of the inflaton. 
Interesting question would be to address a problem where 
there exists a hidden sector which is responsible for super- 
symmetry breaking in that sector, and also responsible 
for mediating supersymmetry breaking gravitationally to 
the observable sector. In such a case the gravitino will 
have an effective mass ~ O(TeV) at a low energy scale. 
So, keeping this in mind we may consider a simple toy 
model with a following superpotential 



1 9 

W = -mi$ 
2 



m 9 [Z+{2- V3)M] 



(37) 



where <!> and Z are inflaton and Polonyi multiplets re- 
spectively. We assume that <f> field is responsible for in- 
flation, so we set mi = 10 13 GeV to produce adequate 
density perturbation, while setting m.2 = 10 11 GeV for 
giving an effective mass to the gravitino around O(TeV). 
An interesting discussion regarding this model has been 
sketched in Rcf. |l2). 

A serious difficulty which immediately arises is that 
one derives a set of coupled equations for the helicity 
±1/2 gravitino and other fermionic degrees of freedom 
Tgl. It has been shown in Ref. that in a global 



supersymmetric limit, this set of equations is reduced 
to a coupled set of equations for the goldstino and the 
transverse combination of the fermions. This suggests 
that there exists a mixing between the goldstino and the 
transverse combination of the fermions. As a result one 
cannot describe the goldstino in a mass eigenstate, and 
thus, it is also difficult to estimate the evolution of their 
number densities. There are many technical difficulties 
because there are essentially two time scales in the prob- 
lem. The first one is related to the fact that the effec- 
tive mass scale of the bosons oscillating and exciting the 
fermionic modes, and, the other one is related to the 
mixing between the goldstino and the transverse combi- 
nation of the fermions (for details, we refer the readers 
Refs. ||[n|). In general one can derive a relationship be- 
tween the two time scales, but this is a non-trivial task 
and we do not have enough tools to address this problem. 

For the above superpotential Eq. (|37|), the inflaton and 
the Polonyi sectors have only gravitational interactions. 
The fermionic components <f> and z have masses mi and 
zero respectively in the global supersymmetric limit. The 
goldstino in this model is a linear combination of the 
fermionic components from both the sectors. As long as 
the energy density is dominated by the inflaton field, the 
helicity ±1/2 gravitinos essentially behave as an inflatino, 
because the mass contribution to the goldstino from the 
Polnyi sector is much smaller ~ O(TeV). This particular 
case is quite interesting and we can analyse the decay of 
the gravitino by assuming that the gravitinos are created 
from the vacuum fluctuations due to the inflaton oscilla- 
tions, whose energy density is dominating the Universe. 
The helicity ±1/2 gravitinos produced during preheating 
will essentially decay because they are essentially the in- 
flatino components and so their couplings are determined 
in the same fashion as that of the inflaton. 

However, the energy density in the inflaton sector is 
decreasing in time, and, when the Hubble expansion ~ 
H < O(TeV), the z component dominates the goldstino. 
Usually, the mixing between the inflatino and z is mini- 
mal and Planck mass suppressed, so, the fermions which 
are produced during preheating will decay again in the 
form of inflatino and cause no trouble for nucleosynth- 
sis, yet there is a finite probability to mix the fermionic 
states and conversion of inflatino to the fermionic part- 
ner of the Polonyi field. Though, we shall not discuss 
this possibility in this paper. One can also imagine that 
the oscillations in the Polonyi sector are also induced at 
H w O(TcV). Once, z field starts oscillating, one might 
expect supersymmetry is broken by the oscillations in 
z direction also, and, as a result gravitinos can as well 
be excited. One may also suspect that the late produc- 
tion of the helicity ±1/2 gravitinos will dominate and 
the problem of gravitino decay still persists. The suspi- 
cion is not fully correct because the number density of 
the helicities ±1/2 and ±3/2 are more or less equal now. 
This is because superpotential contribution to the mass 
of the fermionic component of the Polonyi field is very 
small ~ O(TeV) and the only time-varying scale is due 
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to time- varying mass of the gravitino ~ e zz / 2m2 |VF|/M 2 . 
The presence of the Planck mass suppression prohibits 
explosive production of the gravitinos at late times, so 
especially in the model we have considered, the late time 
production of hclicity ±1/2 can not be very abundant. 
But, now the problem could be much more severe, be- 
cause these gravitinos with both the helicities are pro- 
duced much later, and their effective masses are also very 
small roughly of the order of TeV. This leads to extremely 
slow decay rate of these gravitinos which may cause a 
problem to the Big Bang nucleosynthesis. This picture 
is similar to late production of gravitinos discussed in 
Refs. [^2|. Furthermore, the oscillating Polonyi field leads 
to an even more serious problem, i.e. the moduli prob- 
lem, of which there is no satisfactory way out. 

Finally, we mention and also pointed out in Ref. Jl2| , 
that if the fermionic components mix freely, the inflati- 
nos can be converted to z (which is the field eventu- 
ally eaten by the gravitino). This presumably occurs 
around the time when contributions to supersymmtcry 
breaking from the inflation sector and the Polonyi sec- 
tor become comparable. This problem is analogue to 
the neutrino flavor conversion and the relevant question 
is to ask the conversion probability. As mentioned, we 
beileve that an effcicient conversion will not take place 
for the Polonyi model. An efficient conversion neverthe- 
less results in a large abundance (i.e. comparable to the 
abundances which are produced during preheating) of z 
fcrmion, on top of what is produced due to oscillations of 
the Polonyi field[j]. We notice that if the inflatino decays 
before H w O(TeV), then the abundance of the inflati- 
nos prior to conversion will decrease leading to a smaller 
abundance for z (and consequently helicity ±1/2 graviti- 
nos) even after an efficient conversion. The quantitative 
analysis is beyond the scope of this paper and we leave 
that for future investigation. 

III. CONCLUSION 

Our main result of this paper is to show that the in- 
flatino coupling to the matter field is similar to that of the 
helicity ±1/2 gravitinos. This merely confirms that the 
gravitino interaction with the supercurrcnt actually leads 
to the same interactions as that of the inflatinos when the 
amplitude of the inflaton oscillations is small \<f>\ -C M, 
under the assumption that the helicity ±1/2 component 
of the gravitino behaves as a goldstino for a momentum 
larger than the gravitino mass in a time-varying back- 
ground. Then we have argued that the production of 
helicity ±1/2 states of the gravitino, especially for mod- 
els where supersymmetry breaking scale is dominated by 



^This is the abundance of z fermion which will eventually 
determine the abundance of helicity ±1/2 gravitinos. 



the inflaton energy scale, cannot be considered as a threat 
for nucleosynthesis. Their overproduction can be easily 
understood from the presence of a second derivative of 
the superpotential with respect to the super fields in the 
equation of motion for the hclicity ±1/2 gravitinos. This 
gives rise to an effective mass for the helicity ±1/2 grav- 
itinos, which is equivalent to the mass of the fermionic 
component of the inflaton, known as inflatino. The state- 
ment is true only if the inflaton sector has a single mul- 
tiplet. In some sense helicity ±1/2 states eat the mass 
of the goldstino, which is related to the inflatino by an 
appropriate phase. These states remember their gold- 
stino nature and this is the reason why they are produced 
very efficiently compared to the helicity ±3/2 states. In 
this paper we have argued that the same goldstino na- 
ture come into rescue the late decay of the helicity ±1/2 
gravitino. It has been argued by many authors that the 
helicity ±1/2 gravitinos effectively behave like goldstino 
just after couple of inflaton oscillations. This, together 
with a requirement that the inflaton must decay to give 
a successful nucleosynthesis, leads to an efficient decay 
of the goldstino, or, the helicity ±1/2 gravitinos. Thus, 
they must not survive until nucleosynthesis, and hence 
they should not be considered as a threat to nucleosyn- 
thesis. This argument holds perfectly well for a single 
chiral field where the goldstino is inflatino with some ad- 
ditional phase. However, the extension of this argument 
in some models where there are more than one sectors 
of supersymetry breaking can be made applicable, pro- 
vided supersymmetry breaking scale is still dominated by 
the inflaton energy density. Such a situation can arise if 
there exists a Polonyi field in the hidden sector, which we 
have briefly discussed. However, we still lack a complete 
formal tools to explore all possibilities such as mixing 
between the fermionic components of the inflaton sector 
and the Polonyi sector. This can in principle change the 
abundance of the helicity ±1/2 component of the grav- 
itinos and a detailed study is certainly required in this 
direction. 

It is important to note that the above discussion does 
not apply to the helicity ±3/2 gravitinos. The production 
of these states during preheating is always Planck mass 
suppressed and their existence is also independent of the 
goldstino, so they decay quite late. Due to time varying 
nature of their masses it is always hard to estimate their 
decay rate. It is also true that the helicity ±3/2 states are 
in general produced in less abundance than helicity ±1/2 
states, however, their abundance cannot be neglected as 
pointed out in Refs. |m,0|. For a single multiplet they 
are the only genuine threat to Big Bang nucleosynthesis. 
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